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ABSTRACT 

Rotation is thought to be a major factor in the evolution of massive stars - especially at low 
metallicity - with consequences for their chemical yields, ionizing flux and final fate. Deriving the 
birth spin distribution is of high priority given its importance as a constraint on theories of massive 
star formation and as input for models of stellar populations in the local Universe and at high rcdshift. 
Recently, it has become clear that the majority of massive stars interact with a binary companion 
before they die. We investigate how this affects the distribution of rotation rates, through stellar 
winds, expansion, tides, mass transfer and mergers. 

For this purpose, we simulate a massive binary-star population typical for our Galaxy assuming 
continuous star formation. We find that, because of binary interaction, 20i 10 % of all massive main- 
sequence stars have projected rotational velocities in excess of 200 km s -1 . We evaluate the effect of 
uncertain input distributions and physical processes and conclude that the main uncertainties are the 
mass transfer efficiency and the possible effect of magnetic braking, especially if magnetic fields are 
generated or amplified during mass accretion and stellar mergers. 

The fraction of rapid rotators we derive is similar to that observed. If indeed mass transfer and 
mergers are the main cause for rapid rotation in massive stars, little room remains for rapidly rotating 
stars that are born single. This implies that spin down during star formation is even more efficient 
than previously thought. In addition, this raises questions about the interpretation of the surface 
abundances of rapidly rotating stars as evidence for rotational mixing. Furthermore, our results allow 
for the possibility that all early-type Be stars result from binary interactions and suggest that evidence 
for rotation in explosions, such as long gamma-ray bursts, points to a binary origin. 
Subject headings: Stars: early-type — Stars: massive — Stars: rotation — Binaries: close — Binaries: 
spectroscopic — Galaxy: stellar content 



1. INTRODUCTION 

The origins of the distribution of initial stellar masses 
and initial stellar rotation rates is not yet fu lly under- 
stood (e.g. |Krumk31[Illll ISii^etSJIinil). For the 
latter it is a key issue that stars inherit spin from their 
parental molecular clouds, whose specific angular mo- 
mentum is orders of magnitude l arger than what ca n 
be placed in a rotating star (e.g. rBodcnhci merl Il995l ). 
Based on the observations of low-mass pre-main sequence 
and massive main-sequence stars, it appears that most 
stars reach the zero-age main sequence rotating signifi- 
cantly below their maximum possible rotation sp eed (e.g. 
lHartmann fe Staufferl H981 Huang et all [20Toh . How 
angular momentum is expelled during th e star forma- 
tion process consti tutes a classic problem dMestell 119651: 
McKee fe Ostrikcr 2007; K rumholz et all l2009t iLarsonl 



201(1) . Observed rotation rates of stars on the zero-age 



main sequence provide important constraints for the the- 
ory of star formation. Unfortunately, massive stars be- 
come visible only once core hydrogen burning is already 
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well underway (|Yorkd[T986| ) such that their initial rota- 
tional velocity distribution is not observationally acces- 
sible. 

The embedding of massive stars in the early 
phases of their lives is indeed unfortunate as 
the initial rotation rate is thought to be a fun- 
damental parameter determining their evolutionary 
fate, comparable to their initial mass and metallic- 
ity. Evolutionary models predict that rotation has 
major consequences for the core hydrogen-burnin: 
phase dMaeder fe Mevnetl I2000bt iHeger fe Lan, 



l- burning 
^[2000f 



pnase llMaecter & Mevnetl I^UUUrjt liieger fe Tangerl izuuut 
Hirschi et al.1 120041 lYoon fc Langerl 120051: iBrott et al.l 



2011bt iPotter et al.1 12012 bl : lEkstrom et alll2012[ ). Fur- 
thermore, there is growing evidence that in a fraction of 
massive stars, the final collapse and explosion is governed 
by rapid rotation, giving rise to hyper-ener getic super- 
novae and long-duration gamma-ray bursts (lYoon et al.1 
[20061 iWooslev fc He gen 120061 : IGeorgy et al.H2009D . ~ 

Several processes may affect the stellar rotation rates of 
massive main-sequence stars after they are born, causing 
deviations from the initial values. First, magnetic brak- 
ing may cause them to spin down. In low-mass main- 
sequence stars magnetic breaking is ubiquitous, due to a 
dynamo process operating in their conve ctive envelopes 
(|Skumanichlll972t ISoderblom et al.lll993D . The fraction 
of more massive stars that show evidence for a magnetic 
field appears, however, to be smaller than ab out 15% 
(|Donati fc Landstreetll200ll iWade et al.ll2012aft 

Second, during core hydrogen burning massive stars 



2 



De Mink ct al. 



expand by about a factor of three. It appears, however, 
that a corresponding spin down of the surface layers is 
prevented by an analogous contraction of the stellar core 
and by efficient transp ort of angular momentum from the 
core to the envelope (jEkstrom et al.l l2008bt iBrott etahl 
l2011aj ). independent of the detailed treatment of the in- 
ternal angular momentum transport processes. The ef- 
fects of main-sequence expansion on the surface rotation 
rates thus seem modest as well. 

Third, even non- magnetic massive stars lose angular 
momentum due to their winds ()Langerlll998l ) . The winds 
of massive main-sequence stars, which are fairly well un - 
derstood (jKudritzki fc Pulsl 120001 iMokiem et all [20071) . 
are s tronger at higher mass and metal licitv (iVink et al 



2001) and for more rapid rotation ([Friend k, Abbott 



19861 ). While for the majority of massive stars this wind 
induced angular momentum loss can be neglected — for 
example, stars of LMC composition below 14 M Q lose less 
than 10% of their angular momentum this way ([Langerl 
120121) — it can be important for O-stars and rapid rota- 
tors. 

Finally, — likely most important for the population of 
massive stars as a whole — are changes in the stellar spin 
due to close binary interaction. In recent decades, it has 



formed in close binary systems ( 


Mason et al.lll998l I2009t 


ISana & Evansll2011l iSana et al.l 


2012[). In these systems. 



the spin rate of bot h components can be drastically af- 
fected by tides ('e.g. lZahnlll975l:IHutlll98ltlde Mink et all 
2009al). mass transfe r (IPacketl 119811: IPols et all 119911: 



Petrovic et all 120051: iDervisoglu et all 120101 e.g.) , or 



when the two stars mer ge fe.g.lPodsiadlowski et allll992l : 
iTvlenda et aIll20Tl . ISana et all (12012ft showed that 
about two-thirds of all massive main-sequence O-stars 
are expected to strongly interact with a companion. In 
the majority of cases, after such interaction, either the 
binary is not expected to be recognized as such due to a 
large luminosity ratio and a large period, or there is no 
binary any more due to merging or the break-up of the 
bina ry as consequence of the su pernova explosion of one 
star (jde Mink et al.ll201lLl2012ft . Consequently, any mas- 
sive star population that is not extremely young contains 
a considerable number of apparently single stars whose 
spins are strongly altered by previous binary evolution. 

In this paper we focus on the effects of binarity on 
the evolution of the rotation rates of a population of 
massive main-sequence stars. We do so by construct- 
ing binary population synthesis models that account for 
the latter three effects: the changing moment of iner- 
tia, wind losses, and — most importantly — binary in- 
teraction. It is our aim, after presenting our method 
and assessing the uncertainties, to compare our results 
with observed distributions of ro tational velocities of 
massive mai n-sequ ence stars (e.g. IMokiem et"ai~l 12006 : 
Daflon et al.ll2007l: IWoIff et al.ll2008HHunter et all l2008a: 



Penny fc G ies 2009t iHuang et alil2010L P. L. Dufton et 



al. 2012, submitted, O. Ramirez-Agudelo et al. 2012, in 
preparation). This will provide the first quantitative as- 
sessment of the difference between the initial rotational 
velocity distribution and the observed, present-day dis- 
tributions. 



2. CODE 



We employ a rap i d bin a ry ev olutionary code binary_c 
(jlzzard et all 120041 120061 I2009D tha t makes use of fit- 
ting formulae jHurlev et al.l |2000|) to stellar models 
(jPols et alJll998[ ) to describe the structure of single stars 
as they evolve as a function of their mass, age and metal- 
licity. To study the effects of binary interaction by tides 
and winds and through mass transfer and mergers we 
use the evolutionary al gorithms originally d eveloped by 
iTout et all (fl997t) and iHurlev eTaLl ([200l . These al- 
gorithms approximate the properties and evolution of a 
star after mass loss or mass accretion by switching to the 
prescriptions of a star with the appropriate mass, core 
mass and relative age. These approximations enable us 
to follow the evolution of a particular binary system up 
to the remnant stage in less than a CPU second, which 
is needed to compute extensive grids of models that span 
the multi-dimensional parameter space characteristic for 
binary systems. In addition it allo ws us to explore th e 
effects of uncertain parameters (cf. Ilzzard "elani2009l) . 
For the purpose of this study we updated and extended 
various aspects of this code, which are described in the 
Appendix. A brief overview is given in the remainder of 
this section. 

Throughout this paper, we present results in terms of 
the 'rotational velocity' with which we refer to the tan- 
gential velocity due to rotation at the equator of the star, 
w ro t ■ We account for the flattening of a rotating star as- 
suming Roche geometry. In this approximation, the 
ratio of the radius at the equator and the polar radius 
becomes as large as R cq /R p = 3/2 when the star ap- 
proaches the Keplerian limit, which is given by 



K = JGM/Rf 



oq,K' 



(1) 



where G refers to the gravitational constant, M is the 
stellar mass and i? q,K is the equatorial radius of the 
star that rotates at the Keplerian limit. Equivalently, 
the Keplerian rotational velocity is given by 

vk = K -R cq ,K (2) 

To model deformation by rotation we use the fact that 
the polar radius is hardly affected by rotation. Even 
when the star rotates at the Keplerian limit the polar 
radius deviates by less than 2% from its non-rotating 
value for stars with masses between 3-20 Mq and it de- 
viates b y less than 5% for sta rs with masses between 1- 
60 Mq |Ekstrom et al.ll2008bl ). This property allows us 
to use the stellar radii i?* in our code, which are based on 
non-rotating stellar models, as a good approximation for 
the polar radii of rotating stars, i? P (w) ~ -R P (0) = i?*, 
where we use the shorthand notation w = Q/Q^. The 
equatorial radii can then be obtained by computing the 
shape of the Rochclobe as a function of the rotation rate, 
as is further explained in the appendix. To account for 
gravity darkening in stars rotating near their Keplerian 
limit, we assume that the maximum rotational velocity 
that can be observed is limited to a fraction /nark of 
the K eplerian rotation rate. Based on iTownsend et all 
(|2004l ) we adopt /dark = 0.7. Although this treatment 
is very simple, i t repro duces the detailed simulations of 
ITownsend et al.l (|2004l ) accurately enough for our pur- 
pose. Varying this parameter affects the extend of the 
high velocity tail in our distribution of rotation rates, but 
it does not have a significant on our main predictions. 
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Processes affecting all stars — We account for 
mass loss via stellar win ds usin g prescr iptions by 
iNieuwenhuiizen fc de Jagerl (1199C ) and iVink et all 
(|2001h as described in iBrott et all (|2011a| ). We treat 
the enhancement of stella r winds due to rotation as in 
IMaeder fc Mevnetl (|2000aD . As justified bv lGeorgy et al.l 
( 20111 ). we assume mass loss through stellar winds to 
be spherical when computing the associated angular 
momentum loss. Details can be found in the appendix. 
Our models do not account for effects of magnetic 
fields on angular mom entum loss by stellar winds (e.g. 
lUd-DouiaT ct al. 2009|). We discuss the possible effects in 
Section [5] 

As the star evolves along the main sequence, the outer 
layers expand while the core contracts. We account for 
changes in the mo ment of ine r tia us ing fitting formulae 
based on models bv lPols et al.l (|1998[) . see appendix. The 
internal rotational profile is approximated assuming rigid 
rotation, which is a reasonable approx imation for main- 
sequence stars fe.g. IBrott et al.ll2011al ). 



Processes affecting stars in close binaries — We include the 
effect of tides on the stellar sp ins and the stellar orbit 
(|Zahnlll977t IHurlev et aLl I2002D and the transfer of an- 
gular momentum during mass transfer via an accretion 
disk or direct impact of the accretion strea m onto the 
surface (jUlrich fc Burgerlll976t lPacket|[l98l[) as detailed 
in the Appendix. We assume that the initial stellar spins 
are aligned with the orbit. One of the main uncertainties 
in the treatme nt of binary interact ion is the efficiency of 
mass transfer (|de Mink et al.|[2007l e.g). In particular, it 
is unclear how much stars can accrete after reaching the 
Keplerian l imit. In ou r standard simulatio n s we follow 
iPaczvnskil ()1991l ) and iPopham fc Naravanl (|1991[ ) who 
argue that the accretion disk regulates the mass and an- 
gular momentum flux through viscous coupling allowing 
the star to continue to accrete. We do however limit the 
accretion rate by the thermal rat e of the accret i ng sta r as 
described in the Appendix and in lHurlev et al.l (|2002[ ). In 
section l5.2.2l we discuss the effect of this assumption. To 
account for systems that come in contact during the mass 
transfer phase we consider a critical mass ratio, g C rit, as 
motivated in the Appendix. We assume that the stars 
come in contact when M acc /Md on < <Zcrit, where Md on is 
the mass of the Roche lobe filling star and M acc is the 
mass of the companion. In our standard simulations we 
adopt (/crit.MS = 0-65 if the donor star is a main-sequence 
star, Qcrit.HG = 0.4 when the donor fills it Roche lobe 
when crossing t he Hertz sprung ga p. For more evolved 
donor we follow lHurlev et al.l fl200 2|). Th e effects of these 
assumptions are discussed in Sect. l5.2Tl We assume that 
a merger product efficiently loses the excess angular mo- 
mentum, such that it rotates at near Keplerian rotation 
when it has settled to its thermal equi librium structure 
(jSills et al.ll200l iGlebbeek et al.l [20091 ) . To account for 
mixing of fresh hydrogen into the c entral regions for ac - 
creting stars and mergers we follow lHurlev et al.l (I2002T) . 
We assume that the convective core size adapts to the 
new mass and we describe this process using improved 
prescriptions for the effe ctive mass of the convective core 
(jGlebbeek fc Pols! 120081) . see Appendix. To be able to 
investigate the effect of the uncertainties related stellar 
mergers we implement a parametrized prescription that 
allows us to explore extreme assumptions. In our stan- 



dard simulations we assume that a fraction [Ai oss = 0.1 of 
the total system mass is lost during the merger and that 
a small fraction fi m i K = 0.1 of the envelope is mixed into 
the convective core, see Appendix. The effect of these 
assumptions is discussed in Section [5j 

2.1. Simulating a stellar population 

To investigate the effect of binary interaction on the 
distribution of rotation rates of early-type stars, we sim- 
ulate a stellar population including binaries under the 
assumption of continuous star formation. We approxi- 
mate the distribution of the initial masses of single stars 
and the primary stars of binary systems by 



/ Ml (Mi)dMi ocMf^dMi, 



(3) 



where Mi denotes the ma s s of the prima ry star and a = 
2.35 ± 0.7 (jSalpeterl [19551 : IKroupall200l"l) . Even though 
one can argue whether an initial mass function that was 
derived for single stars can be applied to binary systems 
we note that the effect of this assumption on our results 
is small, see Section [5] 

For the initial distribution of mass ratios, orbital peri- 
ods and the c orresponding bina ry fraction we use the dis- 
tributions bv lSana et al.l (I2012D based on extensive mon- 
itoring campaigns of O stars in nearby yo u ng clusters 
(iDe Becker et all [200l iHillwig et all [200l ISana et al.l 
120081 120091 IRauw et al.l 120091: ISana et al.1 1201 ID . After 
correcting for biases, the authors find a binary fraction 
of 69 ± 9% for systems with mass ratios between 0.1 and 
1 and orbital periods between 10 015 and 10 3 5 days. The 
distribution of orbital periods is described by 



f p {P) dlogPoc (logPf dlogP, 



(4) 



where it = — 0.55 ± +0.2. The distribution of mass ratios 
is described by 

f q {q)dq(xq K dq, (5) 

where q = jM\ , M<i denotes the mass of the secondary 
star and k = —0.1 ± 0.6. Because these distributions 
are derived for a very young population, we assume that 
they well approximate the initial distribution functions 
of binary parameters. In our standard simulations we 
adopt a = 2.35, tt = —0.5 and k = 0. The effects of 
varying these parameters is discussed in Section [5j see 
also Table [U 

We do not include the effect of binaries with compan- 
ions outsi de the range of orb ital periods and mass ratios 
covered by Sana et al. (2012) because the fraction of such 
systems is poorly constrained. As a result we may under- 
estimate the effect of binaries. Systems with mass ratios 
more extreme than 0.1 may produce rapid rotators of 
early type if they interact during the main sequence of 
the primary and merge. 

To reduce the dime nsions of the param eter space we as- 
sume circular orbits. IHurlev et al.l (|2002f) show that tides 
circularize the binary orbit just before the onset of mass 
transfer, implying that the effect of the initial eccentric- 
ity on the further evolution is small. By not considering 
the effect of very wide systems with high eccentricities 
we may underestimate the fraction of stars that interact. 

The initial rotation rate of stars in a binary system 
does significantly affect the evolution of the binary sys- 
tem, except in extreme and rare cases fe.g. lde Mink et al.l 
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l2009al) . The reason is that the stellar spins at birth typ- 
ically contribute at most a few percent to the total an- 
gular momentum budget of the system. 

Since the rotational velocity distribution of massive 
stars at the zero-age main sequence is poorly constrained 
(cf. Sect. 1), we choose a very simple one that allows us 
to clearly demonstrate the effects of binary interaction. 
We adopt a uniform distribution for u rot in the range 
— 200 km s -1 , effectively assuming that stars are born 
with low to moderate rotation rates. While the upper 
limit of 200 km s -1 is not physically motivated, it will 
allow us to unambiguously identify the binary contribu- 
tion to the formation of rapid rotators. Since the initial 
stellar spins are negligible compared to the amount of 
angular momentum exchanged during mass transfer or 
merger, they do not affect the rotation rates of stars af- 
ter these interactions. Whereas in future work, it may be 
considered to adapt the initial spin distribution to repro- 
duce the properties of suitable observed samples, this is 
beyond our present scope. To compute the distribution 
of projected rotation rates fsini, where i is the incli- 
nation angle of the binary system we assume that the 
orientation of the binary orbits are random in space 

Unless an observational campaign is designed to detect 
binaries, many companions stars will remain undetected. 
We assume that only the rotation rate of the brightest 
star is measured in this case. Therefore, when construct- 
ing the simulated distribution of rotation rates we only 
include the most luminous main-sequence star of each 
binary system. 

For our standard simulation we adopt a metallicity of 
Z = 0.008, appropriate for the Large Magellanic Cloud. 
This metallicity is also considered representative for star 
forming regions at a redshift 1-2, at the peak of star 
formation in the Universe. The effect of metallicity is 
discussed in Sect. 14.11 

We derive the distribution of rotation rates for systems 
that are brighter than 10 4 and 1O 5 L0, respectively. To 
put this in perspective, in our models 10 4 Lq corresponds 
to the luminosity of a 8.5-12 Mq main-sequence star, de- 
pending on whether we take the model at zero-age or at 
the end of the main sequence. Similarly, a luminosity cut- 
off of 10 5 Lq corresponds to stars with masses in excess 
of 20-28 Mq. Effectively, the first group is dominated 
by early B-type stars and the second group by O-type 
stars We refrain at this stage from applying other cut-offs 
such as criteria based on temperature, color or spectral 
type, since our predictions for the temperatures are less 
reliable than those for the luminosities. 



3. THE EVOLUTION OF THE ROTATIONAL VELOCITY 
FOR INDIVIDUAL SYSTEMS 

As a star evolves, its rotational velocity is affected by 
various processes, for example as a result of angular mo- 
mentum loss through stellar winds. While a single star 
can only lose angular momentum as a result of mass loss, 
a star in a binary system may either lose or gain angu- 
lar momentum as it interacts with its companion. Even 
when angular momentum is conserved, the rotational ve- 
locity of a star can alter as a result of changes in the 
stellar interior. These effects are discussed in Section |3~T1 
The effect of binary interaction is discussed in Section ^. 21 
and 13.31 
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Figure 1. The evolution of the Kcplcrian rotational velocity for 
stars of different initial masses during their main- sequence evo- 
lution. Labels show the relative age as a fraction of the main- 
sequence lifetimetMgj This diagram is constructed using stellar 

radii bv IPols et a l. (1998) for a metall icity of Z = 0.008, using 
the prescription by [Hurlcv ct al. (200fj). We note that the time- 
dependent radii for stars more massive than about 40 Mq are very 
uncertain. In the most massive stars the decrease of the Keplcrian 
velocity wit h time may well be more severe than shown here. See 
Scction l3.1l for a discussion. 

3.1. Effect of changes in the stellar structure 

Stars expand over the course of their main-sequence 
evolution. Although one might expect intuitively that 
the rotational velocity decreases as the star expands, in 
practice this is not the case. During the main-sequence, 
contraction and consequent spin-up of the core counter- 
act the effect of the moderate expansion of the envelope. 
In this context it is instructive to investigate how the mo- 
ment of inertia, / ~ kR 2 , changes as the star expands. 
Here, k denotes the square of the effective gyration ra- 
dius, which depends on the internal density profile. If we 
approximate k ~ R~^, or equivalcntly 



dink 
dhiR' 



(0) 



we find that the exponent £ varies only slightly during the 
main sequence with typical values of 1.5-2. In the case 
of rigid rotation, we can now express how the rotation 
rate f2, the rotational velocity u rot and the ratio of the 
rotation rate and the Keplerian rate oj / uk change as the 
star expands, 



dlnR * 
din v rot 

dhiR ~ 4 

dinn/n K _ 1 

~~dh[R~~ 2' 



(7) 
(8) 
(9) 



In other words, given typical values of £, we find that 
the decreasing gyration radius compensates for the ex- 
pansion of the star such that the rotation rate Q de- 
creases only slightly as the star evolves (Eq. [7]), since 
d In Cl/d In R « —0.5-0. The rotational velocity at the 
equator, u rot , increases slightly (Eq. [5J. 

Most interestingly, the last expression shows that stars 
naturally evolve towards the Keplerian limit (Eq. O , if 
the amount of angular momentum loss is small and inter- 
nal angular momentum transport between the core and 
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Figure 2. Four annotated examples of the effect of binary interaction on the stellar rotation rates of main-sequence stars for systems with 
different initial mass ratios (left versus right panels) and different initial orbital periods (upper versus lower panels) for Z = 0.008. Each 
panel shows the evolution of the equatorial rotational velocity v TO t of the primary and secondary star (thick lines) as well as the Keplerian 
rotational velocity % (thin lines) as long as the stars are on the main sequence. Shading highlights the phases during which one of the 
stars is rotating more rapidly than 200kms — 1 . See Sect I3.21 for more information. 

envelope is efficient. An important implication of this 
is that when a star reaches the Keplerian rotation rate, 
it remains rotating near the Keplerian limit in the ab- 
sence of an efficient angular momentum l oss mechanism. 
We ref er to the excellent discussion by lEkstrom et alJ 
(|2008b| ) , which describes this effect in detailed models of 
single stars that allow for differential rotation. 

Figure [T] shows the Keplerian rotational velocity «k 
as a function of the initial mass of a star at different 
stages during the main sequence. In zero-age main- 
sequence stars the Keplerian velocity increases monoton- 
ically with initial stellar mass. The Keplerian velocity 
drops as stars evolve and expand. The largest change oc- 
curs during the final stages of main-sequence evolution. 
The change in radius of more massive stars during the 
main sequence is larger, resulting in a more significant 
drop in the Keplerian velocity. As a result the Keple- 
rian rotational velocity at the end of the main sequence 
is around 400 km s^ 1 with only a weak dependence on 
the stellar mass. Note that the projected rotational ve- 
locity for such a star accounting for gravity darkening is 
considerably smaller, by a factor of /dark (sin i) w 0.55. 

We note that the amount by which the massive 
stars expand over the main sequence is not well con- 
strained. The models on which this diagram is based 



(|Pols et al.1 119981 ) have been calibrate d to the radii of 
eclipsing binaries of intermediate mass (|Pols et alJH997t 
iSchroder et al.lll997l ). In these models the expansion of 
massive stars over the course of the main sequence is 
limited to a factor of 2-3. In contrast, the models by 
iBrott et al.1 (|2011a| ). which adopt a larger value for the 
overshooting parameter, predict expansion by a factor 
3-5 for stars in the mass range 5-30 Mq. While the pre- 
dictions based on models with both codes show excellent 
agreement at zero-age, the Keplerian velocities at the end 
of the main sequence are s maller by about a factor of two 
in the IBrott et all (|2011al ) models. 

Furthermore we note that in the most luminous stars 
the effects of radiation pressure can not be ignored and 
the Keplerian limit should be considered as an upper 
limit to the physical maximum. 

3.2. Examples of the spin evolution of binary systems 

In Fig. [2] we depict four examples of the evolution of the 
rotational velocity for main-sequence stars in a binary 
system. In each example we assumed an initial mass for 
the primary star of 20 M Q , a metallicity of Z = 0.008 
and initial rotational velocities of 100 kms -1 for both 
stars. In the top row we show systems with initial orbital 
periods of P = 25 days. These systems are so wide that 
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the primary star fills its Roche lobe only after it leaves 
the main sequence as it expands during its hydrogen shell 
burning phase. In the bottom row we assumed an initial 
orbital period of P = 2.5 days. In these tight systems the 
primary star fills its Roche lobe as a result of expansion 
on the main sequence. 

The panels on the left show examples in which the 
initial mass of the secondary is comparable to that of the 
primary, M2/M1 = 0.75. In these examples one or more 
phases of mass transfer eventually lead to spin-up of the 
companion star. In the panels on the right we adopted a 
more extreme initial mass ratio, M2/M1 = 0.25. In these 
systems the onset of mass transfer brings the stars into 
contact and they merge. Only in the short period case are 
the two stars that merge both still main-sequence stars. 
After the merged star regains its thermal equilibrium it 
is expected to continue to burn hydrogen in the center. 

The phases during which one of the stars is rotating 
more rapidly than 200 km s -1 are highlighted in Figured] 
with grey shading. Notice that during the major part of 
this phase the rapidly rotating star is single or appears 
to be single. The clearest example is shown in panel (d) 
where the rapidly rotating star is the product of a merger 
between the two stars. In panels (a) and (b) the rapidly 
rotating star is the spun up secondary The primary star 
has lost its hydrogen envelope and is hard to detect, as 
a result of its reduced mass, its low luminosity and the 
wide orbit. When the primary star finishes its nuclear 
burning and explodes it is likely to disrupt the system, 
leaving the rapidly rotating secondary behind as a single 
star. 

Besides the drastic effects of mass transfer and coa- 
lescence, Fig. [2] illustrates the effects of other processes 
that have a more subtle effect of the stellar rotation rates, 
which we describe below. 

Changes in the stellar structure — The effect of changes 
on the stellar structure as the star evolves can be ob- 
served, for example, in panel (a) of Figurc[2] During the 
first 9.5 Myr of the evolution of this system both stars 
reside well within their Roche lobe and their evolution 
is similar to that of single stars. The rotational velocity 
of both stars remains roughly constant during this phase 
(cf. Sect. 13. ip . When the primary star approaches the 
end of its main-sequence life time, its expansion acceler- 
ates. This leads to the decrease of the rotational velocity 
that is visible in panel (a) at an age of 8-9,5 Myr. The 
expansion is also responsible for the decrease of the Ke- 
plerian rotational velocity with time. 

Stellar wind — The effect of angular momentum loss by 
stellar winds is small during the major part of the main 
sequence in these examples. The winds become stronger 
toward the end of the main sequence, which together 
with the expansion, contributes to the decrease in the 
rotational velocity of the primary discussed above. How- 
ever, spin down by winds does play a significant role for 
the massive, rapidly rotating stars that can be produced 
as a result of mass transfer. This effect can be seen most 
clearly in panel (c) for ages of 12-15 Myr. The secondary 
quickly spins down reducing its rotational velocity by 
about a factor of three. 

Tides — Tidal interaction tends to synchronize the ro- 
tation of the stars with the orbit in systems where the 



separation between the stars is comparable to the stel- 
lar radii. The systems depicted in the upper and lower 
panels have initial separation in the order of 120 R@ and 
40i?o respectively, slightly less for the systems on the 
right due to the smaller mass ratios. Tides do not play a 
significant role during the main sequence of the stars in 
the upper panels, but they are responsible for the grad- 
ual increase in rotational velocity that can be seen in the 
lower panels for the primary star during the first 8 Myr. 
The orbital period remains nearly constant during this 
phase. The primary star is kept in corotation with the 
orbit as it expands, which implies that the rotational ve- 
locity gradually increases. In this example, the primary 
expands by just over a factor of two before it fills its 
Roche lobe, resulting in a rotational velocity in excess of 
200 kms" 1 . 

In panel (c) of Figure [2] tides arc also important for 
the secondary star, during the first mass transfer phase, 
which starts at about 8 Myr. The secondary star spins 
up as it accretes mass and angular momentum. However, 
the tides quickly spin the star down to synchronous ro- 
tation. This system experiences a mass transfer phase, 
which last almost 4 Myr, during which both stars remain 
in synchronous rotation while the orbit gradually widens. 
Around 12 Myr a second rapid phase of mass transfer sets 
in, i.e. case AB, as the primary star leaves the main se- 
quence and expands during hydrogen shell burning. As a 
result of the high mass transfer rate and the fact that the 
orbit widens, tides are no longer effective in preventing 
the accreting star from spinning up. 

3.3. Effect of the initial separation and mass ratio 

To further illustrate the effect of the initial binary pa- 
rameters, we depict in Figure [3] the evolution of the ro- 
tational velocity of the brightest main-sequence star in a 
binary system as a function of the initial orbital period. 
For this example we adopted an initial rotation rate of 
100 kms -1 , a metallicity of Z = 0.008 and initial masses 
of 20 and 15M Q . The color shading indicates the equa- 
torial rotational velocity of the brightest main-sequence 
star in each system. Initially this is the primary star, but 
after mass transfer the secondary becomes the brightest. 

Mass transfer with a main-sequence donor (Case A) — In 
short period systems the expansion of the primary star 
during its main-sequence evolution is sufficient to make 
it fill its Roche lobe. These systems experience a phase 
of slow mass transfer that can last for several Myr. Mass 
transfer typically occurs via direct impact onto the sur- 
face of the secondary. Tides keep both stars synchro- 
nized, which prevents the secondary from reaching very 
high rotation rates. This phase ends when the two stars 
come into contact and merge (for periods P < 2 days) 
or when the primary star leaves the main sequence (for 
periods P 2-5 days) triggering a new mass transfer 
phase (Case AB). Both cases are expected to lead to the 
formation of a massive rapidly rotating star. In certain 
cases the rapidly rotating star experiences an additional 
spin up phase, when the primary fills its Roche lobe again 
as it expands during He shell burning (Case ABB). This 
effect is visible in Figure [3] around 14 Myr for systems 
with orbital periods of 4-5 days. 

The rapidly rotating stars resulting from such close 
binaries are typically more massive than the original pri- 
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Figure 3. Rotational velocity (color shading) as a function of initial orbital period and time for the brightest main-sequence star in a 
binary system. We adopt initial masses of 20 and 15 Mq, an initial rotational velocities of 100 km s -1 and a mctallicity of Z = 0.008. As 
the stars evolve along the main sequence their rotational velocity is altered by stellar winds, internal evolution, tides and most notably 
mass accretion. The vertical dotted line indicates the maximum separation for which this system interacts by mass transfer. The examples 
shown in panel (a) and (c) of Fig. [2] part of this simulation. 

mary. As a result of their high mass and large luminosity, 
their remaining life time is short compared to the rapid 
rotators resulting from wider systems and they experi- 
ence more efficient spin down by their stellar winds, as 
can be seen in Figure [3] 

The type of evolution described in this section is rep- 
resentative for systems with comparable initial masses 
for the primary and secondary star. In systems with 
more extreme mass ratios, the secondary star is predicted 
to swell up and fill its Roche lobe as well directly after 
the onset of mass transfer. If the two stars evolve into 
a contact situation it is expected that mass is lost via 
the outer Lagrangian point, effectively draining angular 
momentum from the system. This drives the stars into 
deeper contact. We assume that such systems will merge 
to form a rapidly rotating single star. 

Mass transfer with a post-main sequence donor (Case B) 

- Systems wider than about 5 days interact when the 
primary star expands after core hydrogen exhaustion as 
it crosses the Hertzsprung-Russell diagram on its way 
to become a red super giant. This happens during 
hydrogen-shell burning (for periods P w 5 — 100 days) or 
after the ignition of helium (for periods P w 100 — 1000 
daysj3- The period ranges quoted here are only rough 
indications, they depend on the mass of the donor star 
and on metallicity. 

Mass accretion in these systems occurs typically via 
an accretion disk. The companion is quickly spun 
up approaching the Keplcrian velocity after accreting 
only a few percent of it s own mass. According to 
iPopham fc Naravanl (| 19911 ) mass accretion continues af- 
ter the star has reached Keplerian rotation as viscous 
processes transport excess angular momentum outward. 

8 Note that part of thi s group is indicated as Case C in the classi- 
cal classification scheme (jKippciihatm & Wciecrt 1967; Lautcrborn 
119701) , which is based on models for lower mass stars at solar metal- 
licity that reach giant dimensions before helium ignition. In our 
example helium is ignited at smaller radii. Since this group behaves 
similar to Case B systems we discuss them in this section. 



The outer edge of the disk is truncated by tides and feeds 
the excess angular momentum back to the orbit. 

The mass transfer rate is higher for the wider systems, 
in which the donor is more evolved. In our standard sim- 
ulations, the mass accretion rate is limited by the ther- 
mal rate, M*/tkh, of the accreting star, where tkh is the 
Kelvin Hclmholtz timescale. This is independent of the 
separation. As a result we find the general trend that 
tighter systems evolve more conservatively resulting in 
more massive secondaries in comparison with wider sys- 
tems. The more massive secondaries have stronger stel- 
lar winds which cause them to spin down more quickly in 
comparison to less massive secondaries. This effect can 
be seen in the steepness of the color gradient in Figure [3l 
Compare for example the change of rotation rate with 
time for a rapid rotator resulting from a binary with an 
initial orbital period of about 5.5 days with the rapid 
rotators resulting from wider systems. 

Surprisingly, the remaining main-sequence lifetime of 
the secondary is roughly independent of the amount of 
accreted mass or orbital period. This is a consequence 
of two counteracting effects. The higher stellar mass im- 
plies shorter evolutionary time scales, i.e a reduction of 
the remaining life time. At the same time mass accre- 
tion results in mixing of fresh hydrogen into the central 
regions, extending the remaining life time. 

Systems with extreme mass ratios are expected to en- 
ter a phase of common envelope evolution soon after the 
onset of mass transfer. The tighter systems result in a 
merger, but the product is a post main-sequence star, ei- 
ther a blue or red supergiant. Although such mergers are 
very interesting, in this paper we focus on main-sequence 
stars, so we ignore these in the further discussion. In 
wider systems where the binding energy of the envelope 
of the primary is smaller and the momentum and energy 
in the orbit are larger, the envelope can be ejected. This 
process is so fast that it is not expected to significantly 
affect the secondary star. So, we assume its rotation rate 
is unaffected. The only exceptions are systems that are 
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Figure 4. Projected rotation rate distribution for a population of main-sequence stars assuming continuous star formation. The top and 
bottom row show results including only stars or systems brighter than 10 4 and 10 5 Lq . The left panels show the cumulative distribution 
function indicating the fraction of stars with projected rotational velocities larger than 200 km s -1 , F r . The center and right panels show 
the full distribution function and a zoom-in highlighting the contribution of the various binary products. See Sect|4]for further explanation. 



just wide enough to avoid a merger, but in which the 
ejection of the envelope results in significant shrinking of 
the orbit, such that tides lock the orbit of the secondary 
to the orbit of the naked core of the primary star. Such 
systems experience a second phase of mass transfer when 
the primary expands during helium shell burning. How- 
ever, these systems are rare. They do not significantly 
contribute to the population of rapidly rotating stars re- 
sulting from binary evolution. 

Mass transfer with supergiant donor (Case C) — In the 
widest binaries, P w 1500-3000 days in our simulations, 
the donor star develops a convective envelope before it 
fills its Roche lobe. When it loses mass we assume it 
reacts by expanding, r esulting in very high m ass trans- 
fer rates (see however IWoods fc Ivanoval 120111 ) . As the 
Roche lobe typically shrinks, this leads to a run-away 
situation resulting in a common envelope phase. In our 
simulations these systems eject the envelope, leaving the 
secondary star and its spin relatively unaffected. 

A common envelope phase is expected for a very wide 
range of mass ratios. Only systems with initial mass 
ratios close to one can avoid this phase, see for example 
IClaevs et alJ (pflTl . 

4. THE DISTRIBUTION OF PROJECTED ROTATIONAL 
VELOCITIES 

We simulate the distribution of projected rotational 
velocities for a population of main-sequence stars con- 
sisting of single stars and stars in binary systems by 
assuming continuous star formation and a uniform dis- 
tribution of low to intermediate rotational velocities at 



birth w rot < 200 km s" 1 , as described in Scct.[2j The cen- 
tral panels of Figure |4] show the resulting distributions, 
/(usini), using a bin size Awsini = lOkrns -1 . The 
upper and lower panels show our results for brightness 
cut-offs of 10 4 and 1O 5 L , respectively (Sect. HJ. The 
corresponding cumulative distributions, 

v sin i 

F(vsini) = /(usini) Ausini (10) 
o 

arc shown in the panels on the left. 

For a luminosity limit of 10 4 L©, the resulting vsvui- 
distribution is bimodal. The majority of the stars have 
low to intermediate rotational velocities, whereas nearly 
one fifth has rotational velocities in excess of 200 kms' 1 , 
i.e. larger than the maximum of the assumed distribu- 
tion at birth. The first group consists predominantly of 
single stars and members of binary systems that have 
not interacted through Roche-lobe overflow. The rapid 
rotators consist almost exclusively of stars that are the 
product of Roche-lobe overflow. The top right panel of 
Figure 0] shows that this group consists mainly of stars 
that gained mass and angular momentum by mass trans- 
fer. A smaller number are stars merger products, i.e stars 
result from the coalescence of two main-sequence stars in 
a contact binary. 

For stars above 10 5 L©, the shape of the distribution 
changes as can been seen in the two central panels. The 
distribution shown in the upper panel shows two sep- 
arated components. Instead, the rapid rotators shown 
in the lower panel form a tail or plateau that merges 
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with the main component at lower v shxi. This difference 
is the result of angular momentum loss through stellar 
winds, which is stronger for the brighter stars. Even 
though the brightest stars are included in both distribu- 
tions, the distribution in the upper panel is dominated 
by the less bright stars, since the initial mass function 
favors lower mass stars. 

A further difference between the upper and lower plot 
in the central panel of Fig. 4 is that the bright star dis- 
tribution extends to higher rotational velocities. The 
reason is that many binary products rotate near their 
Keplerian limit, which is larger for the brighter, more 
massive stars (cf., Fig. 1). 

Not all binary products have large projected rotational 
velocities. Apart from rapid rotators observed at small 
inclinations, i.e. near pole-on, there is a small contri- 
bution from systems that are undergoing mass transfer, 
indicated with the shorthand "Algols" in Figure 2] This 
group consists almost entirely of short-period systems 
that are undergoing a slow, i.e. nuclear timescale, Case A 
mass transfer (cf. Sect. 3). The rotational velocity re- 
flects the rotation rate of the accreting star, which has 
become the brightest star in the system. Even though 
this star is accreting mass and angular momentum, tides 
arc efficient enough, during this slow phase of mass trans- 
fer, to prevent it from spinning rapidly. These systems 
are easily detected as binaries. As the donor star fills 
its Roche lobe these systems are likely to show eclipses 
or at least ellipsoidal variations. Examples of such sys- 
tems are the 30 semi-detached sy stems in the sample of 
double-lined eclipsing binaries bv lHilditch et~aLl (I2005T) . 

The distribution of the rotational velocities of single 
stars closely resembles the adopted uniform birth distri- 
bution. The apparent bias towards low rotational veloc- 
ities is only partially due to spin down by stellar winds 
and changes in the stellar structure. The main reason is 
the projection effect caused by the random distribution 
of spin axis of the stars. 

The rotational velocities of stars in binary systems that 
have not yet interacted by Roche-lobe overflow can be 
affected by tides. The effect of spin up by tides can be 
seen in the yellow shaded area as a tail of stars rotat- 
ing between 200 and about 300 km s -1 , which is most 
pronounced in the lower central panel. Tides are also 
responsible for the flattening of the yellow curve for pre 
Roche-lobe overflow systems below 100 kms" 1 , as they 
counteract the spin down towards the end of the main se- 
quence which is imposed by stellar winds and stellar ex- 
pansion. Even though tides may also lead to spin down, 
they do not produce very slow rotators, since slow syn- 
chronized rotation implies wide orbits and tides are no 
longer effective in wide binaries. 

We have included mergers separately in Figure HI to 
show how the distribution may change if mergers behave 
differently than we have assumed in our simulations (see 
Sect. 15.21 for further discussion). 

4.1. Metallicity dependence 

The metallicity affects the distribution of rotational 
velocities via its effect on the stellar wind and on the 
stellar structure. In Figure[S]wc show the normalized dis- 
tribution of rotational velocities for main-sequence stars 
brighter than 10 5 Lq for a range of metallicitics. We 
assume here that the initial mass function, the binary 
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Figure 5. The distribution of projected rotational velocities 
/(usini) for different metallicitics ranging from Z = 0.001 (blue 
line) to 0.03 (red line) for systems brighter than L > 10 5 Lq 
containing at least one main-sequence star. The inset gives the 
color code and summarizes the effect of metallicity by showing the 
fraction of stars with project ed r otational velocities larger than 
200 kms -1 , F r , i.e.. See Sect. [47X1 for further information. 

fraction and the binary distribution functions arc not 
dependent on metallicity. 

At the lowest metallicitics shown, the distribution con- 
sists of two distinct components with broad peaks at 
vsini ~ 100 and 400 kms -1 . It resembles the distri- 
bution shown in the top panels of Figure 0] for 10 4 Lq, 
although the low metallicity distribution extends to ve- 
locities as high as 600 kms -1 . At the highest metallicity 
the two components merge to one smooth distribution 
with a strong peak at vsini ~ Okms . 

The inset in Figure [5] shows the fraction of stars with 
projected rotational velocities in excess of 200 kms -1 , 
F r , as a function of metallicity, Z. This fraction is in the 
range 13-22%, with larger values for smaller metallicities, 
and a stronger metallicity dependence for more massive 
stars. These trends are mainly the result of the depen- 
dence of the stellar winds on metallicity and brightness, 
and the corresponding wind-induced angular momentum 
loss. 

Additionally, main-sequence stars tend to be more 
compact at lower metallicity, where hydrogen burning 
via the CNO cycle operates less efficiently. More com- 
pact stars have higher Keplerian rotational velocities for 
the same mass. This is reflected in the extent of the high 
velocity tail in Figure [5] which increases with decreasing 
metallicity. Furthermore, the smaller radii imply some- 
what more compact binary systems, which raises the im- 
portance of tidal interactions at low metallicity. 

5. EFFECTS OF UNCERTAINTIES IN THE ADOPTED 
ASSUMPTIONS 

In this section we take advantage of the speed of our 
rapid binary evolutionary code and population synthesis 
routines, which enables us to evaluate the impact of our 
assumptions. This concerns the initial distribution func- 
tions (Sect. 15. ip and assumptions about the physics of 
binary evolution (Sect. 15. 2\i . We quantify the effects by 
comparing the fraction of stars with projected rotational 
velocities in excess of 200 kms -1 , F r , that results from 
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Table 1 

Impact of the assumed initial binary fraction and t he i nitial 
distribution of binary parameters (as defined in Sect. |27TJ) on the 
fraction F r of stars with vs'mi > 200 km s -1 . For discussion sec 
Sect. 



initial distribution 


variation 


F r (%) 


standard* 




18.8 


orbital period 


> 7T > -1 


13.1 - 24.7 


mass ratio 


-1< K < 1 


14.8 - 21.3 


primary mass 


1.65 < a < 3.05 


17.0 - 21.3 


binary fraction 


0.5 < / bin < 1.0 


14.6 - 23.9 



Note. — *tt = -0.5, k = 0.0, a = 2.35 and / bin = 0.7 



the various assumptions. We use our standard simulation 
for L > 10 4 Lq as reference. 

5.1. Initial distribution functions 

The effect of uncertainties in the adopted initial distri- 
butions is relatively small, each affecting the fraction of 
rapid rotators only by a few percent in absolute terms, 
see Table [1] The initial orbital period distribution is the 
largest cause of uncertainty accounting for a variation 
of ±5.8%. When we compare a distribution that is flat 
in logP, i.c Opik's law (tt = 0) to a distribution that 
strongly favors short period systems (tt = —1), we find 
that F r increases from 13.1 to 24.7%. The latter distribu- 
tion favors systems that interact while the primary star 
is still on the main sequence. In our simulations, these 
systems produce rapid rotators through mass transfer 
and mergers that are on average more massive, and thus 
brighter, than rapid rotators produced in wider binaries. 

Uncertainties in the adopted distribution of mass ra- 
tios are responsible for a variation of just over ±3% 
when we vary the power law exponent from k = — 1 
to k = 1, which gene r ously allows for the uncertainties 
quoted by iSana et al.l ([2012:). A distribution skewed to- 
wards systems with comparable masses for the compo- 
nents (ft = 1) favors systems that produce rapid rotators 
through nearly conservative mass transfer. The spun up 
mass gainers are typically brighter than the original pri- 
mary and they experience significant rejuvenation. Both 
effects result in a larger fraction of rapid rotators in the 
population. A distribution that favors systems with ex- 
treme mass ratios (k = —1) results in a larger number 
of mergers. While mergers in short period systems con- 
tribute to the fraction of rapid rotators, mergers resulting 
from wider systems do not as we consider these to become 
post main-sequence objects, explaining the reduction of 
the fraction of rapid rotators 

The rapid rotators in our simulation result from stars 
that gain mass, apart from the small contribution of 
tidally spun up stars. On average a gain in mass im- 
plies a larger luminosity. Stars that were initially not 
bright enough to be included in our distribution can be- 
come sufficiently bright after accreting. Changing the 
relative number of lower mass systems, by adapting the 
initial mass function, increases the number of potential 
systems that can produce rapid rotators. However, the 
effect is small: varying the steepness of the initial mass 
function of the primary star by a — 2.35 ± 0.7 leads to a 
change of about ±2%. 

The effect of varying the initial binary fraction can be 
obtained by rescaling the results. We account for the 



fact that the fraction of stars born as single stars in the 
population F s is not equal to the fraction of single stars 
at birth. The fraction of rapid rotators F r for any given 
initial binary fraction /t,i n is given by 

Fr(M a ) = F;^( F ^ + F:\^-) \ (11) 

•/bin \ -/bin <f bin / 

where superscript * refers to the values obtained in our 
standard simulation (see Fig. E} and Ft, = 1 — F s . Con- 
sidering a population in which only half or all stars are 
born as binaries with mass ratios and orbital periods as 
specified in Sect 12.11 results in variations in F r of about 
±5%. 

5.2. Uncertain physics 

In this section we address uncertain physical assump- 
tions regarding the treatment of mass transfer, the for- 
mation of contact systems and mergers. 

5.2.1. Do contact systems merge? 

Mass accretion drives the companion star out of ther- 
mal equilibrium causing it to expand when the mass 
transfer timcscalc is short com pared to the t hermal 
timcscale of the companion (e.g. INeo et al.l[l97l This 
can lead to the forma tion of a contact binarv ifPolsll 1 9941 
IWellstein et al.ll2001| ). However, the parameter range for 
which contact occurs is uncertain as it depends, among 
other things, on the poo rly constrained spec ific entropy 
of the accreted material (|Shu &: Lubow|[l98lj) . 

When contact is established, it is not clear whether this 
always leads to a common envelope phase and a merger. 
If the secondary star overfills its Roche-lobe by only a 
small amount, the star may detach again as the star re- 
stores thermal equilibrium ()de Mink et al.ll2007D . How- 
ever, the outcome depends on mass and angular momen- 
tum loss during this phase. If the stars come into deep 
contact, mass loss from the system through the outer 
Lagrangian point efficiently drains angular momentum 
from the system and drives the stars deeper into con- 
tact. Mass loss through stellar winds or bipolar outflows 
can, in principle, have the opposite effect. 

In our standard simulation we account for the forma- 
tion of contact due to rapid mass transfer by considering 
the mass ratio q = M^jMx at the onset of Roche- lobe 
overflow. If the mass ratio is smaller than a critical value 
we assume that contact is established and results in a 
merger. If the mass ratio is larger, we continue to com- 
pute the accretion of mass and angular momentum onto 
the secondary star, and assume that contact is avoided 
as long as the thermal equilibrium radius of the mass 
gainer is smaller than its Roche radius. Below we dis- 
cuss the impact on our results of varying the adopted 
critical mass ratio systems that interact during the main 
sequence q C rit,MS and systems in which the donor star is 
a post main-sequence star crossing the Hcrtzsprung gap 

9crit,HG- 

Contact involving two main-sequence stars — Systems that 
interact during the main sequence eventually produce a 
rapid rotator cither through a merger or through spinup 
by during mass transfer regardless of their initial mass 
ratio. Therefore the fraction of rapid rotators is not very 
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sensitive to the adopted critical mass ratio for merg- 
ing <?crit,MS- In our standard simulation we adopted 
9crit,MS = 0.65. If we vary g CT it,MS in the range 0.8-0.25, 
the fraction of rapid rotators changes by only ±0.5%. 
On the one hand, mergers are typically more massive 
and more luminous than the secondary star would have 
been had the system not merged. As a result, in our 
brightness limited sample the rapidly rotating stars re- 
sulting from mergers originate from binary systems with 
a wide range of initial masses and mass ratios. On the 
other hand, we find that the remaining relative lifetime 
of mergers is typically much smaller than that of rapid 
rotators produced through mass transfer. Mergers con- 
tain the helium rich core of the original primary and even 
though some fresh hydrogen is mixed into the central re- 
gions, their remaining lifetime is limited. In contrast, the 
core of the secondary is still hydrogen rich at the onset 
of mass transfer in comparison to the core of primary. 
In the population, the longer remaining lifetime of the 
secondaries compensates for their lower brightness. 

Contact involving Hertzsprung gap donor — For post main- 
sequence interaction the formation of contact introduces 
a larger uncertainty on the predicted number of rapid 
rotators. If a binary system mergers after the primary 
star has left the main sequence, the product is expected 
to become a post main sequence object, a blue or red 
super giant. Hence, we do not include such stars in our 
distribution of rotation rates. In our standard simula- 
tion we assume that such systems merge if the mass ra- 
tio at the onset of Roche lobe overflow is smaller than 
9crit,HG = 0-4. Detailed binary evolutionary models show 
that the parameter range for th e formation of contac t 
may be substantially larger (e.g. iWellstein et all 120011 ). 
For the extreme assumption that all these systems merge 
we find an absolute change of —5.4% for the fraction of 
rapid rota tors. To obtain this e stimate, we mimic the 
results by IWellstein et al.1 ()2001f ) by assuming that all 
systems with initial orbital periods larger than 30 days 
and mass ratios smaller than 0.65 merge. Even though 
this assumption affects a very large fraction of the initial 
parameter space, the effect is limited because it mostly 
affects systems that would produce relatively low mass 
rapid rotators. Taking the opposite extreme assumption 
that none of the Case B systems merge (g C rit,HG = 0.0) 
does have a significant effect on the fraction of rapid ro- 
tators, for the same reasons. 

5.2.2. Can stars continue to accrete after they are spun up 
to Keplerian rotation? 

IPacketl (|1981f ) argued that a star fed by an accre- 
tion disk reaches Keplerian rotation after accreting just 
a few percent of its initial mass. It remains unclear 
whether a star rotating at this s peed can cont i nue to 
acc rete mass. In the mo dels by iPetrovic et al.l (|2005f ) 
and Ide Mink et all (|2009b|) it was assumed that accretion 
ceases when the star reaches the Keplerian limit. Only 
in short period systems tides can spin the accreting star 
down such that it can undo ubtedly resume accreting a 
substantial amount of mass. IPopham fc NaravanI (|1991f ) 
argue that viscous coupling between the star and the 
disk can govern an inward flow of mass and an outward 
flow of angular momentum such that also mass gainers 
in wider binaries may continue to accrete while rotat- 



ing near or at the Keplcrian limit. While we are not in 
a position to resolve this debate, we have the tools to 
investigate the impact of both extreme assumptions on 
the distribution of rotation rates. In our sta ndard simu- 
lation we followe | Popham fc NaravanI ()1991h . Following 
IPetrovic et al.l (|2005l ) instead does not reduce the num- 
ber of stars that are spun up, but it does reduce the mass 
and therefore brightness of these stars. As a result we 
find that the fraction of rapid rotators in our brightness 
limited sample changes by —11%. 

5.2.3. Mixing and mass loss during mergers? 

Mergers constitute the least understood phase of bi- 
nary evolution. In particular, the amount of mass loss 
and the amount of mixing is uncertain. In our standard 
simulations we assume that a fraction fi\ oss = 0.1 of the 
total system mass is lost during a merger and that a frac- 
tion /x m i x = 0.1 of the envelope is mixed into the core. 
Changing the amount of mass loss affects the number 
of mergers that are bright enough to be included in the 
sample, but the effect is small. Varying fi\ oss in the range 
0-0.25 results in a change that is not larger than ±0.4%. 
Changing the amount of mixing primarily changes the 
remaining lifetime of the merger product. Varying /Lt m ; x 
in the range 0-1 to simulate the effect of no additional 
mixing versus complete mixing of the core and the enve- 
lope we find a change of —0.8% and ±2.5%, respectively. 

5.2.4. Magnetic fields? 

Recent searches have revealed l arge scale magnetic 
fields in several O type stars (e.g. I Wade etldl I2012bf ) 
with strengths large enoug h to affect their spin down 
time ([Ud-Doula et al.1120091) . Several recent findings may 
be relevant for our study. The apparent single O star 
HP 148937 ha s shown to host a large-scale magnetic field 
(IHubrig et al.ll2008l MTU. INaze et al.l f2010t MaAeetal\ 
l2012cD , and is a stron g candidate for being a merger 
product (|Langer|[2012f ). It is surrounded by a bipolar 
nebula that may have been ejected during the merger 
event. Furthermore, the recent detection of a large 
scale magnetic field in the mass gainer of Plaskett's star 
(|Grunhut et al.l [2012ft implies that such fields can exist 
in a star shortly after an accretion event, and may even 
indicate that such fields can be generated as its result. 
In this context it is interesting to mention the tentative 
det ection of a magnetic fi eld in the O s tar of Cyg-Xl 
by iKaritskava et ail (|2010l ) see however Bagnul oet al.l 
(12012J), and for HD 153919 (IHubrig et al.H201ll) . whichTs 
the O star co mpanion to the ~ 2.5 M & compact object 
4U1700-37 (cf. lClark et alll200l . 

If the production of a large-scale magnetic field would 
be a commo n by-product of a m a ss transfer or merger 
event (e.g. iFerrario et al.l 120091 : iTutukov fc Fedoroval 
I2010D or if could be gen erated as a result of a rotation- 
ally driven dynamo (e.g. IPotter et alj|2012al) . we would 
have unde restimated the spin down times of such prod- 
ucts (e.g. iMevnet et al.l I2011D . Since this scenario is 
speculative at the moment, we refrain from quantifying 
such an effect, but it might reduce the number of rapid 
rotating stars produced by binary interaction. It would 
not, however, decre ase the number of bin ary products in 
a given population ()de Mink et al.ll2012D . 

6. DISCUSSION AND IMPLICATIONS 
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Table 2 

A compilation of observational studies listing the percentages of stars, F, with projected rotational velocities larger than 200, 300, and 

400 fans -1 . 



Reference Environment Sp.Type Sample size F >200 kms -i F >300 kms-i ^oo km S -i 
(%) (%) (%) 



Observed samples 



Conti & Ebbets (19771 


MW 





205 


26 


10 





Howarth et al. (19977 


MW 


OB 


373 


9 


4 


1 


Abt et al. (2002) 


MW 


BO-3 


357 


16 


4 


0.3 


Martavan et al. (2006) 


LMC (NGC 2004) 


B 


121 


15 


7 





Martavan et al. (2007) 


SMC (NGC 330) 


B 


198 


27 


9 


1.0 


Hunter et al. (2008a) 


LMC 


OB 


204 


16 


2 





5? 


SMC 


OB 


204 


22 


6 





Pennv & Gies (2009) 


MW 


OB 


97 


22 


10 


1 


LMC 


OB 


106 


8 


9 







SMC 


OB 


55 


13 


5 


0.2 


Huang et al. (2010) 


MW (Clusters) 


BO-9 


695 


28 


7 


0.1 


MW (Field) 


BO-9 


483 


19 


5 





LT-FLAMES Tarantula Survey (VFTS) 












Dufton et al. (2012, subm.) 


LMC (30 Dor) 


09.5-B3 


337 


50 


15 


3 


Rarmrcz-Agudelo et al. (in prep.) 


LMC (30 Dor) 


O 


178 


30 


16 


5 



Predictions from this work 

all 19 11 2 

mergers 5 3 0.5 

other 14 8 1.5 



Based on our simulations we expect 20^ % of massive 
main-sequence stars to rotate rapidly as a result of bi- 
nary interaction. In this section wc discuss implications 
for rotationally induced mixing in stars (Sect. 16. ip . and 
for our understanding of the origin of Be stars (Scct. l6T2|) . 
We place our results in context of the observed distribu- 
tion of rotation rates (Sect. 16. 3p . and we discuss how to 
proceed from here to derive the true birth spin distribu- 
tion of massive stars (Sect. I6.4[) . 

6.1. Implications for rotationally induced mixing 

Rotation has been argued to be a major factor in- 
fluencing the evolution of massive stars. In particular, 
it is thought to induce mixing in the radiative layers 
of m assive stars (jMaeder fc Mevnetl I2000bt iHeger et al.l 
120001 and references therein) . The presence of hydrogen- 
burning products at the surface of O and B type 
stars has generally been interpreted as a signature of 
rotational mixing dGies k Lambertl 119921: IHunter et al.l 
I2008bl: [Pr^billa et al.ll2010T ). IHunter et al.l ^OOSbF arT 
alyzed about one hundred early B-type main-sequence 
stars in the LMC with projected rotational velocities of 
up to ~300 kins -1 . The rapid rotators with enhanced ni- 
trogen abundances in this sample have been considered 
to provide the most direct evidence for rotational mixing 
(jMaeder et al.H2009t IBrott et al.H2011al) . 

Our results question the validity of this interpretation; 
they suggest that a significant fraction, or perhaps even 
close to all of the rap id rotators are th e products of close 
binary evolution (see lAcke et al.ll2008l for a potential ex- 
ception). Detailed models show that mass transfer in a 
binary results in sur face nitrogen enrichment that cover s 
the observed range (jWellstein et al.|[2001l ; lLangerl [2012h . 
If the observed nitrogen enrichment of the rapid rota- 
tors is primarily the result of mass t ransfer, there must 
be less room for r otational mixing (|Langer et al.l [20081 : 
IBrott et al.ll2011bD . 

Additional evidence for the importance of close binary 



evolution in the sample of IHunter et al.l (2008b) come s 
from the star-by-star analysis of iKohler etaiT (12011) . 
They conclude that 10 of the rapidly rotating stars in 
the sample can not be explained by rotating single stel- 
lar models; the observed nitrogen abundance for these 
stars is much lower than expected for their age and their 
projected rotation rate. To reconcile these objects with 
the theory of rotational mixing requires to assume that 
these stars were slow rotators for most of their lives, 
and that they have been spun-up only recently due to 
close binary interaction by non-conservative mass trans- 
fer. We conclude that it is (so far) difficult to single out 
the effect of rotational mixing in samples of massive main 
sequence stars. It turns out that the strategy to circum- 
vent this problem by removing identified binaries from 
the observed sample is not e ffective; it may even ac hieve 
the opposite. As argued in Ide Mink et al.l (|2011[ ) , the 
products of binary evolution typically are or appear to 
be single stars. As most of the stars which are found to be 
a member of a binary system did not yet gain any mass 
from their companion, removing stars with evidence for 
binarity from the sample may in fact increase the rela- 
tive contaminati on of the sample wi th post-interaction 
binary products (jde Mink et al.ll2012l ). 

We conclude that the role of rotational mixing re- 
mains ambiguous at present. This is unsatisfying given 
the large implications of rotational mixing on the evo- 
lution of the progenitors of long g amma-ray bursts 
and pair-instability supernovae (e.g. lYoon elaTl l2006t 
IWooslev fc Heeen 120061 ILaneer et al.l 120071) as well as 
on the chemical and radiative feedback o f massive stars, 
in particular in the ear ly Universe (e.g. Ekstro m et al.l 
l2008aHYoon et al.|[2012l) . Further studies of the relation 
between rotation rates and nitrogen surface abundances 
of large and well-defined samples of stars are essential. 
A first step may be to establish whether the types of 
behavior seen in B stars are also present in O stars. 
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6.2. Evidence for spin-up from Be/X-ray binaries and 
the binary origin of Be stars 

Compelling evidence for spin-up of massive main- 
sequence stars as a result of binary interaction comes 
from the Be/X-ray binarie s, the most fre quent type of 
high mass X-ray b inaries dLiu et al l [20061) . This class, 
first recognized by iMaraschi et al.l ( 19761 ). consists of a 
neutron star in an eccentric orbit around a rapidly ro- 
tating B-main sequence star. The formation of these 
systems can be understood by assuming that the pro- 
genitor of the neutron star transferred mass and angu- 
lar momentu m to its companion before it explo ded as 
a supernova (|Rappaport fc van den Heuvell Il982i) . The 
Be/X-ray binaries constitute the subset of systems that 
remained bound after the explosion. The majority of the 
systems is, however, expected to be disrupted as a result 
of the kick of the neutron star, thereby producin g single 
Be stars (e.g. !Blaauwlll96ll: lEldrklge et al.ll20Tl . 

The idea that a significant fraction, or even all, of 
the Be st ars result from b inar y interaction has been ex- 
lored bv lPols et all ()1991l ) andlvan Bever fc Vanbeverenl 
1997f l. Both studies conclude that binaries provide a sig- 
nificant fraction of single Be stars, but they also find that 
the number of mass transfer systems is not high enough 
to explain all Be stars. However, our models show that 
binaries produce significantly more rapid rotators than 
quoted in these studies. This is partially due to the fact 
that we consider stellar mergers as a channel to form 
rapid rotators. Furthermore, the initial binary parame- 
ter distri butions ad opted in these original s tudies were 
based on|Ab! (|1983f ) and lAbt fc Levvld 19781). We adop t 



distributions based on recent work (jSana et all I20121) . 
which made clear that earlier studies undere stimate the 
frequ ency of short period binaries for O-stars (jSana et al.1 
12012) 

We refrain from predicting the fraction of Be stars pro- 
duced by our models because the precise conditions for 
the Be-phenomenon remain somewhat elusive. It is how- 
ever worth to mention that the total fraction of mergers 
and mass gainers in our standard simulation is 24.1%, 
(cf. the top right panel of Fig. 4). For comparison, the 
fraction of Be stars among early-type, non-supergiant B- 
type stars is 20-30% after correcting for selection effects 
(jZorec fc BriotJ[T997T) . 

6.3. Comparison with observed samples 

A direct comparison between the various observational 
results, summarized in Table [3J and our predictions is 
difficult, due to different selection and methodological 
effects, and potentially different age distributions of the 
respective stellar populations. For example, some ob- 
servational studies attempt to identify and remove spec- 
troscopic binaries while others did not. Furthermore, 
some samples are biased against rapid rotators, for ex- 
ample because emission line stars were not included or 
because of an inherited bias in samples that are based 
on archival spectra. A comparison with our predictions 
should thus also be taken with care. Overall, however, 
we find remarkable trends which appear to be present in 
all samples. 

The fraction of stars with usini > 200kms _1 in the 
observed samples compiled in Table [2] ranges between 
10 and 50%, while our calculations predict numbers of 



the order of 20%, all produced from binaries. Our the- 
oretical number is based on the assumption of a con- 
stant star formation rate. We expect that for certain 
star formation histories, e.g. a starburst, and certain ob- 
servational biases, e.g. in samples focusing on the cluster 
turn-off, the rates may be significantly different over the 
situa t ion captured by our models (e.g. iPols fc Marinusl 
[llMlvanBever fc VanbevererifllM IChen fc Hanll2008L 
Schneider et al. 2012, in prep). 

Considering higher cut-off values for v sin i appears to 
leave less and less room for single stars, even consid- 
ering the error margin in our prediction (cf., Sect. 5). 
Above 300 kms" 1 , only the Dufton et al. (2012, subm.) 
and Ramirez- Agudelo et al. (2013, in prep.) samples, 
both resulting from the VLT-FLAMES Tarantula Sur- 
vey (VFTS), find more stars than we predict. The lat- 
ter works however have removed detected spectroscopic 
binaries from their sample, which account for almost a 
third of their total sample. As can be seen from Fig. 4, 
spectroscopic binaries are expected to mostly contribute 
to the v sin i distribution below 300 kms^ 1 . If we correct 
our predictions for this we obtain -P^ooians- 1 = 16%, 
which is remarkably close to what is observed in the 
VFTS samples. 

This trend continues when considering stars above 
400kms _1 . While the number of observed stars is small, 
the associated percentages from the VFTS samples are 
close to our results. We note that the VFTS samples do 
not suffer from an explicit bias against the most rapid 
rotators, in contrast to the other samples. Whether this 
is sufficient to explain the observational differences with 
the VFTS results or whether there is a genuine difference 
between the rotational properties of stars in 30 Dor and 
stars in less extreme environments would require a care- 
ful assessment of the selection effects and observational 
biases of these other samples. 

The above analysis suggests a decreasing fraction of 
genuine single stars in stellar samples with increas- 
ing v sin i cut-off. It appears possible that there is a 
v sin i threshold above which binary products dominate 
the population. To establish such a threshold, and its er- 
ror bar, requires tailored binary population models which 
reproduce specific observed samples; this is beyond the 
scope of this paper. Nevertheless, from the preliminary 
comparison of our models with the existing data in Ta- 
ble [5] it seems conceivable that this threshold value - 
which will depend on metallicity, star formation history 
and sample biases — could generally be smaller than 
300 km s" 1 . 

6.4. Towards deriving the initial distribution of rotation 

rates 

Deriving the initial distribution of rotation rates is a 
high priority given its im portance as a con straint for 
star-formation theories (e.g. lRosen eTaT]|20l"2l ) and as in- 
put condition for modeling populati ons of massive stars 
both nearby and at high redshift dBrott et all 1201 lbt 
iLevesque etaLll2012t lEldridge fc Stanwavll2009f In the 
literature, it is often assumed that the distribution of ro- 
tation rates of early-type stars closely reflects their birth 
spin distribution, arguing that the rotational velocity of 
massive stars is not expected to change significantly dur- 
ing their main-sequence evolution — which we do confirm 
when considering only genuine single stars. However, we 
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have shown that this assumption docs not generally hold, 
due to the large close binary fraction of massive stars 
(|Sana et al.ll2012t) and the consequences of close binary 
evolution. 

For star clusters or a star-formation rate which is 
strongly peaked in time, the effects may be smaller (e.g. 
for extremely young star clusters) or larger (looking at 
the main-sequence turn-off of star clusters, where all 
stars are close to ending core hydrogen burning). Since 
star clusters that are much younger than the lifetime of 
their most massive stars are difficult to find, it may re- 
main a challenge to derive the true birth distribution of 
the spins of massive stars directly from observations. 

There may be two promising ways to proceed. One 
may try to observe the birth spin distribution of mas- 
sive stars directly in young star clusters, excluding the 
most massive stars. If such kind of studies are success- 
ful, the v sin i-distribution can be compared with those 
of stars in the same mass range, for older clusters, where 
binary effects had time to operate. This could be com- 
pared with results from our method, with tailored input 
functions appropriate for the observed samples, which 
would allow to draw conclusions about the effects of bi- 



nary evolution. Alternatively, our method may be used 
in the future to directly derive the initial spin distribution 
of an observed population, adopting an iterative process 
where this distribution is varied until the observations 
are appropriately reproduced. 
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APPENDIX 
RAPID BINARY EVOLUTIONARY CODE 

In Section [2] we briefly described the rapid binary evolutionary code employed in this study. Here, we provide a 
more complete description with references and document routines that have not been documented elsewhere. 

Stellar wind mass loss 

To accoun t for the radiatively driven winds of massive stars we implement the mass-loss prescription of IVink et al.1 
(|2000l[200lT ). This recipe accounts for the fast increase of the mass- loss rate for stars with temperatures below 22,000 K, 
which is related to t he recombination of Fe IV to Fe III and is commonly referred to as the bi-stability jump. Following 
iBrott et al.l (|2011aD . we perform a linear interpolation in the mass loss rate within 2,500K of the jump, to ensure a 
continuous transition. 

To accommodate the strong mass-loss i ncrease when approaching t h e iHumphrevs fe D avidsonl (|1994f ) limit, we 
switch to the empirical mass loss rate of INieuwenhuiizen fe de Jaeer (119901) at any temperature lower than the 



critical temperature for th e bi-st ability jump, when the IVink et al.l ()2001l ) rate becomes smaller than that from 
iNieuwenhui j zen" fe de Jagerl (fl990t ). This' ensures a smooth transition between the two mass loss prescriptions and 
naturally accounts for the increased mass loss at the theoretical ly predicted second bi-stability jump at ~ 12,500K. 

We adop t Wolf- Rayet mass - loss rates of lHamann et al.l (|1995l ) for naked helium stars and stars with thin hydrogen 
envelopes (|Hurlev et al.ll2000l) . To account for the extreme mass loss of s t ars be yond the Humphreys-Davidson limit 
limit an LBV- like mass loss rate is included as describe d inlHurlev et al.l (120001). M ass-l oss rates during the evolve d 
sta ges of interme diate and low- mass stars arc based on lKudritzki fe Reimersl 1 19781 ) and IVassiliadis fe Wood! (|1993| ). 
see iHurlev et~aT1 (puML 

The prescription bv lVink et all (|2001[ ) accounts the metallicity dependence of radiatively driven winds with a scaling 
factor of (Z/0.02) - 85 , where Z is the metallicity expressed as i nitial mass fraction of elements heavier then helium. 
We apply the same scaling for the INieuwenhuiizen fe de Jagerl (|1990( ) rates. Other rates are assumed not to depend 
significantly on metallicity (e.g. Ivan Loom I2006D . 

Effects of rotation 

Throughout the evolution of each binary system we compute how the stellar spins change as a result of their internal 
evolution and angular momentum loss and gain by tides, mass exchange and stellar winds. 

Internal evolution and the moment of inertia 

As the star evolves its mo ment of inertia , I = k MR changes, where k is the radius of gyration squared. The original 
implementation described in lHurlev et al.l (|2000D assumed constant k = 0.1 during the main sequence. However, even 
over th e course of the main sequence the evolution of the internal structure significantly changes (e.g. lEkstrom et al.l 
l2008b| ). Therefore, we adopt fitting formula (Pols, priv. communica tion) for the evolution and mass dependence of 
the gyration radius based on evolutionary models bv lPols et al.1 ()1998l ). For zero-age main-sequence stars the gyration 
radius squared, ko, is given by 

k ~ c + min{0.21, max{0.09 - 0.27 log M, 0.037 + 0.033 log M}}, (Al) 
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where the correction factor c = except for stars with logM > 1.3, for which c — — 0.055(logM — 1.3) 2 . When the 
stars evolve along the main sequence, the star becomes more centrally condensed and the gyration radius decreases. 
This can be described in terms of the radius R and the radius at zero- age Rq, 



R\ c „^fR^- 01 



k ~ (k Q - 0.025) I — J + 0.025 ^— J , (A2) 
where 

f -2.5 logA/<0, 
C = { -2.5 + 5 log M for < log M < 0.2, (A3) 
{ -1.5 0.2 < log M. 

Although our fits are based on models at a metallicity Z = 0.02, we account for the main effect of metallicity 
through the dependence of the radii on metallicity, as a result of the choice to express k in terms of R/Rq. This 
approximation is sufficient for our purposes. Analytic approximations of the evolution of the gyration radius at more 
evolved evolutionar y stages are included as well, but the differences with respect to the original implementation by 
IHurlev et al.l (|2000f ) do not significantly affect our study. 

We assume that the internal rotational profile can be approximated by the assumption of rigid rotation as a result of 
efficient intern al transport of angular momentum. For m ain-sequence stars this is supported by the findings of detailed 
stellar models ([Ekstrom et a l. 2008b; Br ott et"aLll2011al ) independent of the detailed treatment of the internal angular 
momentum transport processes. 

Deformation by rotation and the Keplerian limit 

Stars become oblate when they rotate near their Keplerian limit, i.e. the rotation rate for which material at the 
equator is no longer bound to the star because the outward centrifugal acceleration balances the inward gravitational 
acceleration. Deviati ons from spherical symmetry have been o bserved directly by interfcromctric studies, for example 
Altair and Alchernar ((Peterson et alJ|2006HCarciofi et a l. 2008). We account for deformation of rotation as described in 
Sect. [21 To obtain the equatorial radiii we compute the shape of the equipotential surface in the Roche approximation. 
Under the assumption th at the polar radius is not affected by rotation this becomes an algebraic equation of the third 
degree (e.g. iMa cder 2009j). We find that the numerical solution obtained by the Newton method can be approximated 
by the following analytic approximation 

log 10 ^(cj)~w (/itan[1.422w] + /2sin[w]), (A4) 
Hp 

where fx = 1.7539 x 10~ 2 and f% = 4.511 x 10~ 2 , which can be evaluated in the timely fashion required for population 
synthesis. We note that there is also some ambiguity in the literature with respect to these two possible definitions 
of 'the fraction of break-up' expre ssed in terms of the the rotation rate, and in terms of the rotational velocity v cq . 
Note that Q/Qk 7^ v cq /v cq ^K (e.g. IMacder 200^). For slow and intermediate rotation rates the fraction of break-up in 
terms of the rotation rate and the fraction of break-up in terms of the rotational velocity deviate by a factor 0.67 since 

J^l = J^-n/n K . (as) 

The interplay of rotation and the stellar wind 

To ac count for the effect of rotation on the mass loss by stellar winds we follow the approach by iMaeder fc Mevnetl 
(|2000al Eq. 4.30). They d erive the latitudi nal mass flux from a star deformed by rotation taking into account the 
effect of gravity darkening (jvon Zeine]lH92l . 

M(u>) ( 1 - r 



M( ) ~ \i-r-f{u)j (A6) 

wher e T is the Eddington factor and a the force multiplier, for which we use the empirical relation by ([Lamers et al.l 
119951 ) valid for logT eff = 3.9 - 4.7, where a varies from 0.15-0.52. Outside this range we assume that a is independent 
of the effective temperature. The dependence on the rotation rotation rates enters through 



/(w) = — 0.198 «-^H . (A7) 
2-rrGpm \ R P J 

The approximation on the right hand side is accurate up to w < 0.8. We compute the Eddington factor, 

r = (as) 

by approximating the opacity k with the electron scattering opacity, k ~ 0.2(1 + X), where X is the mass fraction of 
hydrogen at the surface. Since we do not explicitly follow the surface composition we adopt X = 0.74 for stars that 
still have a hydrogen envelope and X = for naked helium stars. 
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In this prescription the effect of rotation on the mass loss rate is small ([Maeder fc Mevnetll2000al) . Only for T > 0.64 
and logT c ff < 4.3 is a large increase of the mass loss rate predicted. In all cases we limit the stellar wind mass loss 
rate to the thermal mass loss rate Mkh = M/tkr, where the we approximate the Kelvin Helmholtz timescale tkh as 

T KH = 10 7 yr|^, (A9) 

where M, R and L refer to the current mass, radius and luminosity expressed in solar units. 

The radiatively-driven stellar wi nd of a star deforme d by rotation becomes asphcrical as a result of the changes of 
the effective gravity with latitude. iGeorgy et al.l ([2011D computed the net effect of the latitudinal dependence on the 
loss of angular momentum. They find that the deviation from the spherical cases is very small. In the interest of 
computational speed we can therefore safely assume that the specific angular momentum of mass lost through the 
stellar wind is equal to j = 2/3i?pf2. 

Mass and angular momentum loss of stars rotating close to Keplerian rotation 

During their life stars ma y approach the Kepler ian rotation rate as a result of mass accretion iPac ket (1981) or as a 
result of internal evolution. lEkstrom et al.l (|2008 b) . In this case the star must dispose of its excess angular momentum 
at a rate dictated by the process that drives the star to its Keplerian limit. If angular momentum loss by the stellar 
wind or tides is not efficient enough, it is likely that a near-Keplerian disk is formed when the thermal motion of gas 
particles at the equator becomes suffi cient to overcome the escape speed, such that they are launched in orbit around 
the star ([Okazaki k, Neguer uela 200 lj). Viscous coupling of the d isk and the star may help extract angular momentum 
efficiently lPaczvnskilll991t iSills et~al]|2005l: iKrticka et al.ll2011l) . while losing only very little mass, depending on the 
viscosity and the extent of the disk. The disk may be truncated by the torque of the companion star. A further 
possibility is the illumination of the disk which may cause it to flare and disperse. The interplay between these 
processes will determine how much mass the rapidly rotating star can accrete. A detailed treatment of outflowing 
disks within a binary system is beyond the purpose of this study, but we discuss the effects of extreme assumptions in 
Section O 

Treatment of binary interaction 
Mass transfer rate 

We model mass transfer as described in iHurlev et al.l (|2002f ). but with the following adaptations. In the original 
version of the code the rate for stable mass transfer is assumed to be a steep function of the amount by which the star 
overfills its Roche lobe (|Hurlev et al.ll2002l eq. 58). This prescription can introduce numerical instabilities causing the 
star to oscillate in and out of its Roche lobe. Therefore, we replace this prescription by an adaptive algorithm. When 
the star fills its Roche lobe we determine the mass transfer rate by removing a small amount of mass and computing 
the resulting change in radius until the stellar radius becomes s maller than the Ro che lobe radius within a certain 
threshold. The rate is capped by the thermal rate as defined in IHurlev et al.l (|2002l eq. 60 ). As a result, the mass 
transfer rate is governed by the evolutionary changes in the stellar radius and changes in the Roche lobe. 

When a star accretes mass on timescales that are short co mpared to the thermal time scale the accreting star may 
be driven out of thermal equilibrium resulting in expansion (jBensonl 119701 ) . How severe this effect is is uncertain as 
it depends on the poorly constrained specific entropy of the accreted mat erial. Shocks during the accretion process 
may result in the accretion of lower entropy material (jShu &: Lubowlll981l ). Most detailed binary evolutionary codes 
use a simple boundary conditions where material is assumed to be accreted with the specific entropy of the surface of 
the accretor. This leads to severe expansion in binaries with extreme mass ratios implying that a large fraction of the 
systems will come in contact with the likely result of a stellar merger. Due to the practical limitations of the design 
of our rapid code we cannot follow the radii of stars during the accretion process. To account for systems that come 
in contact as a result of this effect we adopt a critical mass ratio for main-sequence binaries, <7 cr it.MSi and assume 
that the stars come in contact when M acc / M& on < q C rit,MS) where Md on is the mass of the Roche lobe filling star and 
M acc is the mass of the companion. This critical mass ratio depends only weakly on the separation (Ide Mink et al.l 
120071 ). We therefore consider two free parameters <7 CT it,MS and g C rit,HG which set the critical mass ratio for systems 
that interact on the main sequence and systems in which the donor star fills its Roche lobe while it is crossing the 
Hertz sprung gap. Our standard value q CV i t ms = 0-65 is based on detailed binary evolutionary models (Ide Mink et al.l 
I2007T) and for 

<7crit,HG we adopt a value of 0.4. For the advanced evolutionary stages we follow IHurlev et al.1 (|2002l ). 

Rejuvenation and the convective core 

If the accreting star is a main-sequence star we assume that it adapts its interior structure to its new mass. This 
implies that the convective core grows, such that fresh hydrogen is mixed to the central regions effectively rejuven ating 
the accreting star. We treat this by adapting the effective fractional main-sequence life time as in lTout et alJ (|T997l ).The 
original implementation implicitly assumed that the size of the convective core is a fixed fraction of 10% the stellar 
mass. This strongly underestimates the core mass for massive main-sequence stars found in detailed models and 
ignores the fa ct that mass fract i on of t he convective core increases with stellar mass. To improve this adopt fit formula 
provided by iGlebbeek fc Polsl (|2008l) for the effective core mass the end of the main sequence as a function of the 
initial stellar mass. We do not account for the shrinkage of the convective core over the course of the main-sequence 
evolution. 
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Transfer of angular momentum accretion 

When mass is transferred through the interior Lagrange point to the secondary star is may either form an accretion 
disk around the secondary star or it may directly i mpact on the surface of th e secondary star. Direct evidence for this 
is found in the Ha profiles in Algol type binaries (jRichards fe Albrightlll999l ) 

To distinguish between direct impact and the formation of an accreti on disk we estimate the minimum i? m i n distance 
between the stream and the accreting star using an analytic fit by (|Ulrich fc Burger] Il976f ) against calculations of 
(|Lubow fc Shulll975l ) to determine whether the stream will hit the star (Pols, priv. communication). iLubow fc Shul 
(|1975[) estimate the specific angular momentum of the impact stream to be equal to \/ GM1.7R min . If the impact 
parameter is larger then the radius of the accretor, the stream is assumed to collide with itself after which viscous 
process lead to the formation of a Keplerian accretion disk. The star is assumed to accrete from the inner radius of 
the disk, where gas particles accrete with the specific angular momentum of a Keplerian orbit with a radius equal to 
the stellar radius y/GMR cq . 

iPacketl (|1981[ ) pointed out that, under these assumptions, the accreting star is very efficiently driven towards the 
Keplerian limit after accreting only a few percent of its ow n ma ss in the case of accretion v ia a disk. In our standard 
models set we will follow the argument of iPaczvns ki (1991) and Popha m fc Naravanl (|199lD who argue that a star ro- 
tating at the Keplerian limit star can keep accreting mass from an accretion disk without accreting angular momentum 
as a result of viscous coupling. 

While the star rotates near the Keplerian limit the net accr eted angular mom entum is slowly reduced below the 
specific angular momentum of a Keplerian orbit as discussed by iColpi et al.1 (|1991| ) . There is some debate whether an 
equili brium situation may be reached at a sub Keplerian rotational velocity, for example in the presence of a magnetic 
field (jGhosh fc Lamblll979l : IPriugk; 1989). If this equilibrium is reached at a very small fraction of the Keplerian rate 
it would affect the results in this study, similar as discussed in Sec l5.2.4l 

Wind accretion 

We consider wind accretion (jBondi fc Hovlel [l944T> with adopting an efficiency factor cibh = 1-5 as detailed in 
I Hurley et al.l (|2002f) . However, in contrast to Hurley et al. we assume that t he specific angular momentum from 
material gained by wind accretion is small and can be ignored (cf. IRuffertl[l99l . 

Treatment of stellar mergers 

A lar ge fraction of massive stars are foun d in binaries that are expected to evolve into contact systems and eventually 
merge |Bensonl H970: Wcllstei n~eTaI1[200l . The mergers that are of main interest for this study are those that result 
from the merg e r of tw o main-sequence stars. We use the extensive grid of detailed binary models evolutionary models 
Ide Mink et al.l ()2007l ) to calibrate the critical mass ratio for close b inaries leadin g to co ntact systems. 

For the treatment of mergers involving an evolved star we follow IHurlev et al.l (f200l . The original prescription for 
mergers involving two main-sequence stars assumes that the merger product is completely mixed and that no mass is 
lost during the process. These assumptions in combination with the original prescription that severely underestimated 
the mass of the convective core for massive main-sequence stars implied that the mergers experience severe rejuvenation 
placing them near the zero-age main sequence. 

We adapt the treatment of mergers between two main-sequence stars with respect to the original implementation to 
better reflect current insights and detailed models. We assume that a certain fraction /ii oss of the combined mass is 
lost from the system when two stars merge. In our standar d simulation we adopt a fraction of 10% in our standard 
simulations in agreement with lLombardi et al.l (|1995l [T996() . Note, however, that these simulations are for direct 
collisions, not for merging binaries. We therefore consider ^ti oss one of the uncertain parameters. We assume that the 
mass lost from the system originates from the stellar envelopes and is not enriched in helium. We assume that, after 
short phase lasting a thermal timescale the merger product will settle and can be described as a stars rotating near 
the Keplerian rate, numerical reasons. 

Whether merger products experience significant mixing is uncertain (e.g. lGaburov et al.l l2008 : ) To mimic the result of 
SPH simulations which show remarkably little mixing occurring when two stars merge we envision in our standard model 
that the core of the most evolved star settles in the center of the merger surrounded by the core of its companion. The 
most evolved star is not necessarily the primary star. We then ass ume that the merger pr oduct develops a convective 
core of a mass that can be approximated by the prescriptions by [Glebbeek fc Pols! ( 2008[ ). Depending on how much 
mass is lost from the system, the new convective core mass may be larger or smaller than the combined mass of both 
original cores. 

We consider the possibility that an additional amount of mass just above the convective core can be mixed into the 
central regions. To this extent we introduce a second parameter ^ m i x whic h expresses the mass of this region as a 
fraction of the envelope mass. We adopt /i m j X = 0.1 as a standard option (Gaburo v et al.l l2008T l . Setting /i m i X = 1 
allows for complete mixing of the merger product after mass loss. We work out the new relative age by considering 
which regions are mixed and following the approach by (jHurlev et al.l l20 02f) to use a simple linear map between the 
core composition and the relative age. 

In this new approach, which uses a more realistic approximation for the core size, the mergers are the relative age of 
the merger product is only slightly smaller than the relative age of the original primary star. In other words we now 
provide a more conservative estimate of the remaining lifetime of the merger products. 
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Common envelope evolution 

To treat common envelope situations we follow IHurlev et al.1 (|200l we use choose a common e nvelope ejection 
efficiency parameter agE = 0-2 in agreement with studies of binary systems in planetary nebulae (|Zorotovic et al.l 
l2010tlDavis et al.ll2010t iDe Marco et alj|201lh . For the parameter describing the binding energy of the envelope Ace 
we use fits to (|Dewi fc Taurisll200(f " and we adopt Ai on i zat i on = 0.5. We note that in the original version of the code 
massive stars are not treated as convective giants when they become red super giants after the ignition of helium. We 
corrected this using the base of the giant branch as a transition point. During common envelope evolution we assume 
that the companion does not accrete mass nor angular momentum. Our results are not sensitive to these assumptions. 
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